Theory of c-axis Josephson tunneling in d x 2 y 2-wave superconductors 



Kazumi Maki and Stephan Haas 
Department of Physics and Astronomy, University of Southern California, Los Angeles, CA 90089-0484 

(February 2, 2008) 

The temperature and angular dependence of the c-axis Josephson current and the superfluid 
density in layered d x 2_ y 2-wave superconductors are studied within the framework of an extended 
Ambegaokar-Baratoff formalism. In particular, the effects of angle-dependent tunneling matrix el- 
ements and Andreev scattering at grain boundaries are taken into account. These lead to strong 
corrections of the low-temperature behavior of the plasma frequency and the Josephson current. 
Recent c-axis measurements on the cuprate high-temperature superconductors HgBa2CaCui+^ and 
Bi2Sr2CaCu20g+i can therefore be interpreted to be consistent with a d x 2_ y 2-wave order parameter. 
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Interfaces such as grain boundaries in the high- 
temperature (high-T c ) superconductors are known to 
have important theoretical, technical, and practical im- 
plications. [1] In a set of early experiments it was shown 
that the critical current between single-crystal surfaces 
depends sensitively on their relative orientation. [2,3] 
More precisely, when two such samples are matched along 
their crystal a-b planes, but twisted with respect to each 
other about the direction of their common c-axis, it was 
found that the critical current decreases rapidly as a func- 
tion of the twist angle a. This strong angular dependence 
has since been observed frequently at grain boundaries in 
the high-T c cuprates. It is a consequence of the d x 2_ y 2- 
wave symmetry of the superconducting order parameter 
in these compounds. [1] 
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FIG. 1. Schematic representation of two single-crystal sam- 
ples, matched along their a-b planes and twisted about the 
c-axis by an angle a. 

In the recent literature, there have been theoretical as 
well as experimental puzzles regarding the precise na- 
ture of c-axis Josephson tunneling in the cuprate super- 
conductors, particularly in I^Si^CaC^Os+i. On the 
theoretical front, several scenarios have been considered 
to model the angular behavior of the Josephson current. 
[4-6] However, in these frameworks it appears difficult 
to reproduce the strong dependence on the twist angle a 
which has been observed in experiment. On the exper- 
imental side, there is a recent controversy arising from 
a c-axis Josephson tunneling experiment in which two 
single crystals were twisted about the c-axis and subse- 



quently joined together. [7] In contrast to earlier exper- 
iments [2,3], the authors found a critical current that is 
very small but independent of a, implying that the su- 
perconducting order parameter in I^S^CaC^Os+s has 
a dominant s-wave component. [7] However, it appears 
that the tunneling area in this particular measurement 
may have been too large to guarantee homogeneity of 
the contact areas. Furthermore, the small magnitude of 
the observed Josephson current suggests a large contri- 
bution from the surface area, which in turn would not 
adequately reflect the symmetry of the bulk supercon- 
ducting order parameter. [8] 

We will analyze a set of more recent experiments, 
which measure the c-axis Josephson current in super- 
conducting Bi 2 Sr 2 CaCu208+5. [9,10] Two single-crystal 
whiskers of this compound with very small widths $ ~ 
10/iTO were joined along their a-b planes. By varying 
their relative orientations the angular dependence of the 
Josephson current along the c-axis was detected. A clear 
fourfold symmetry with respect to the twist angle a was 
observed, consistent with a d x 2_ y 2-wave superconduct- 
ing order parameter for this compound. Surprisingly, 
however, the angular dependence turned out to be much 
more pronounced than one would have expected from 
previous theoretical treatments, based on the Ginzburg- 
Landau model and on the Ambegaokar-Baratoff formal- 
ism for superconductor-insulator-superconductor (S-I-S) 
junctions, which predict a simple cos (2a) dependence for 
d x 2_ y 2-wave superconductors. [11] 

Here we propose an extension of the Ambegaokar- 
Baratoff treatment for the superfluid density and the 
Josephson current along the c-axis of layered unconven- 
tional superconductors, taking into account (i) the coher- 
ent component of the tunneling current, (ii) the d x 2_ y 2- 
wave symmetry of the order parameter, and (iii) the an- 
gular dependence of the transfer matrix element along 
the crystal c-direction. 

Let us first discuss the issue of coherence. In contrast 
to the case of s-wave superconductivity, Josephson cou- 
pling between two unconventional superconductors with 
nodes in their gap functions is possible only if there is 
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a coherent component to the tunneling current. Hence, 
purely incoherent tunneling will not result in a Joseph- 
son coupling between two d x 2_ y 2-wave superconductors, 
as already noted in Refs. [4] and [12], and for the layered 
k-(ET) 2 salts in Refs. [13] and [14]. Coherent tunneling 
between adjacent layers with planar quasiparticle wave 
vectors k and k' can be defined by the condition k \\ k' . 
This is what Ambegaokar and Baratoff called "specular 
transmission" . [11,15] On the other hand, if the tunneling 
is incoherent k ^ k', i.e. k and k' are completely uncor- 
rected. For the sake of the following argument, let us 
consider only these two extreme cases. In real materials, 
coherent and incoherent processes are typically mixed, 
in which case only the coherent component is relevant. 
The Josephson coupling is proportional to (A(</>)A(</>')), 
where (• • •) represents angular averages over <f> and <f/, 
which are the angles between the wave vectors k and k! 
and the crystal a-axis. Hence, when k and k! are uncorre- 
cted (AO) AO')) = (AO)) (AO')) = 0. This argument 
applies not only to d x 2_ y 2-wave systems, but also to other 
unconventional superconductors. In contrast, there will 
be Josephson coupling between two s-wave superconduc- 
tors even when the tunneling is completely incoherent, 
and therefore this is one of the important differences be- 
tween conventional and unconventional superconductiv- 
ity. 

When the transport along the c-axis is due to coher- 
ent tunneling, the temperature dependence of the out- 
of-planc superfluid density p sc is different from that of 
the in-plane superfluid density p s \\. Recent calculations 
performed for anisotropic three-dimensional models have 
predicted a linear decrease with increasing temperature, 
both for p sc and p s ». [16,17] However, within the tunnel- 
ing model, the out-of-plane superfluid density for d x 2_ y 2- 
wave superconductors is given by 
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where t = T/T c is the reduced temperature and A(t) is 
the amplitude of the superconducting gap function. Thus 
the tunneling model predicts a ^-dependence at low 
temperatures, consistent with the out-of-plane superfluid 
density observed in YBa 2 Cu 3 07_,5 [18], Tl 2 Ba 2 Cu0 6 
[19], and k - (ET) 2 Cu(N(CN) 2 )Br [13]. 

More recently, the out-of-plane superfluid density in 
Bi 2 Sr 2 CaCu 2 0s+5 was measured, using a sensitive mi- 
crowave resonance technique. [20] Noting that the square 
of the Josephson plasma frequency is proportional to 
p sc (t), it was observed that the low-temperature behav- 
ior of the c-axis superfluid density is proportional to 
t 5 ~ t 6 . A comparable high power- law has also been 
reported in HgBa 2 CaCui + a. [21] In order to describe 
this behavior, Xiang and Wheatley [22] introduced an 



angle-dependent tunneling matrix element proportional 
to cos 2 (20), which vanishes along the nodal directions. 
This behavior is common to those high-T c cuprates which 
do not contain Cu-0 chains. It is caused by the overlap 
of Cu 4s and O 2p orbitals. [20,22] Introducing an angle- 
dependent tunneling matrix element into the formalism 
for the superfluid density leads to 



Psc(t) 



i5AQ) r' 2 

8A J 



> sin 



, f A(t)sm<p 
> tanh 



2T 

6 



(2) 



. T \ T 
1 - 3.606 — - 443.35 — 

a ; V^o 



While the leading power appears to be t 3 , comparison 
of the coefficients of the t 3 - and i 6 -terms shows the t 6 - 
power to be dominant at all experimentally relevant low 
temperatures. 
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FIG. 2. Out-of-plane superfluid density as a function of 
temperature. The solid line includes the angular depen- 
dence of the tunneling matrix element along the c-direction, 
whereas the dashed line does not. The inset shows the 
low-temperature region. 

In Fig. 2 the out-of-plane superfluid density of a 
d x 2_ y 2-wave superconductor is shown as a function of the 
reduced temperature. Comparing the behavior of p sc (t) 
with and without inclusion of the angular dependence 
in the tunneling matrix element, the difference appears 
minor to the eye. However, by zooming into the low- 
temperature regime (inset of Fig. 2) one observes that 
the extra cos 2 (2</>)-term leads to a strongly reduced tem- 
perature dependence in this regime. A similar difference 
in the low-temperature behavior has been reported in 
measurements of the plasma resonance of various high- 
Tc cuprate compounds (Fig. 3 of Ref. [20]). 

Let us now consider the angular and temperature 
dependence of the Josephson current, as probed us- 
ing Bi 2 Sr 2 CaCu 2 0s+5 cross-whiskers twisted about their 
common c-axis. [10] By applying the Ambegaokar- 
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Baratoff formalism, and assuming that there is a com- 
ponent of coherent tunneling along the c-direction, one 
obtains 
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where w„ = 2nT(n + 1/2) is the Matsubara frequency, 
and t± is the tunneling matrix element along the c- 
direction. The angular dependence of the gap function is 
contained in the factors 

/>=cos(20), /< -cos (2(0 + a)), (4) 

if | cos (20)| > | cos (2(0 + a)) |, 

/>=cos(2(0 + a)), /<=cos(20), (5) 

if | cos (20)| < | cos (2(0 + a)) | . 

The factor in the denominator (1 + (A /t_i_) 2 (|/> | — 
,/<|) 2 ) includes the effect of Andreev reflections [23] at 
the grain boundary. If this factor is neglected, one ob- 
tains the usual Ambegaokar-Baratoff expression, with a 
dependence of the Joscphson current on the twist angle 
very close to cos (2a). Furthermore, one notes that this 
correction of the Josephson current due to Andreev scat- 
tering vanishes if a = 0, i.e. when the two whiskers are 
perfectly aligned. 

Taking into account the additional angular dependence 
of the tunneling matrix element in the above derivation 
leads to a slightly modified expression, analogous to the 
discussion of the supcrfluid density, 
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Let us first discuss the angular dependence of the zero- 
temperature limit of Eqs. 3 and 6, given by 
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where K(k) is the complete elliptic integral of the first 
kind. Here F(0) = 1 for the case of angle-independent 
tunneling along the c-axis (Eq. 3), and F((f>) = |/</>| 2 
when this angular dependence is included (Eq. 6). 

Fig. 3 shows the angular dependence of the c-axis 
Josephson current for various values of the ratio of the 
gap amplitude to the tunneling matrix element Ao/t±. 
Data from the recent cross-whisker measurements on 



Bi2Sr2CaCu20s+<5 by Takano et al. [10] are also shown 
as filled circles. In Fig. 3(a) the tunneling element is as- 
sumed to be constant. Therefore the angular dependence 
is less pronounced than in Fig. 3(b) where an additional 
cos 2 (20) modulation of t± has been included. Compar- 
ing with the experimental data, both models appear to 
fit the strong angular dependence of the measured cur- 
rent reasonably well with A /t± ~ 10 for F(4>) = 1, and 
A /i_L ~ 20 for F(0) = |/</>| 2 . However, the rapid 
decay of I C {T = 0, a) appears to be better captured by 
the model which includes the additional angular modula- 
tion of the tunneling matrix element. In this comparison 
with the experiment we have assumed that the data point 
at a = 0.4887T which violates the otherwise monotonic 
behavior in I C (T = 0, en) is of extrinsic origin, perhaps 
related to inhomogencitics in the contact surfaces. 
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FIG. 3. Dependence of the zero-temperature c-axis Joseph- 
son current on the twist angle a. Solid lines: calculated 
I c iT = 0,a) according to Eq. 7, (a) with F(0) = 1, and 
(b) with F((f>) = |/</>j 2 . From top to bottom the ratio of 
the gap amplitude to the c-axis tunneling matrix element is 
given by A /t± = 0, 1, 2, 5, 10, 20, and 50. Circles: data 
from Ref. 10. 

The most important conclusion from the above con- 
siderations is that Andreev scattering has a crucial role 
in determining the strong angular dependence of the 
Josephson coupling across grain boundaries. As observed 
in Fig. 3, the additional factor (1 + (A /t±) 2 (|/>| - 
|/<|) 2 ) -1 leads to substantial deviations from the simple 
proportionality I c (T,a) oc cos (2a), which one obtains 
without taking Andreev processes into account. This 
observation is also related to the question of whether 
these grain boundaries are S-N-S or S-I-S junctions. The 
present study indicates that they seem to be something 
in between. 

Finally, let us use the above formalism to make some 
predictions regarding the temperature dependence of the 
c-axis Josephson tunneling current for various twist an- 
gles. From the discussion of the angular dependence it 
is expected that the current exhibits peaks at a = and 
7r/2, and vanishes at a = 7r/4 and a = 37r/4. In anal- 
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ogy to the data presented in Ref. [10], we show in Fig. 4 
the temperature dependence for selected twist angles in 
the region a G [7r/4, 7r/2]. The ratio A /ij_ was fixed 
according to the best fits obtained in Figs. 3(a) and (b). 
One observes in Fig. 4 that while the main temperature 
dependence is governed by the prefactor A 2 (t) ~ (1 — t 3 ) 
in Eqs. 3 and 6, the most noticeable differences between 
the two cases, F(<p) — 1 vs. F{4>) = |/</>| 2 , occur in the 
low-temperature regions. Moreover, they are more pro- 
nounced at large twist angles close to a = w/2. Also, 
exactly at a = ir/2 the expression for the Josephson 
current reduces to p sc (t). Comparison of future exper- 
imental data with these curves will aid in distinguishing 
between the proposed models. 
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FIG. 4. Temperature dependence of the c-axis Joseph- 
son current. From top to bottom the twist angle is 
a = tt/2, 0.45tt, 0.4tt, 0.35tt, and 0.3tt. (a) F{4>) = 1 and A/t± 
= 10. (b) F(4>) = |/</>| 2 and A/t ± = 20. 

In conclusion, we have studied the angular and tem- 
perature dependence of the superfluid density and the 
Josephson current across c-axis grain boundaries, such as 
the recently synthesized twisted cross-whisker junctions 
of Bi2Sr 2 CaCu208+,5. For both quantities, p sc {t) and 
I C (T, a), the inclusion of an additional cos 2 (2(f)) angular 
modulation in the tunneling matrix clement is impor- 
tant for accurate modeling of the low-temperature prop- 
erties of nodal cuprate superconductors without Cu-0 
chains. The observed strong angular dependence of the 
Josephson current, deviating substantially from a sim- 
ple cos (2a) behavior, can be accounted for by Andreev 
scattering processes at the grain boundaries. Therefore 
the measured Josephson current across whisker surfaces 
in Bi 2 Sr 2 CaCu208+5 appears to be fully consistent with 
a d x 2_ y 2-wave symmetry of the superconducting order 
parameter. 
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